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Correlations among the optical properties of 
cirrus-cloud particles: Microphysical interpretation 

J. Reichardt, 1 ’ 2 S. Reichardt, 1 ’ 2 M. Hess, 3 and T. J. McGee 2 

Abstract. Cirrus measurements obtained with a ground-based polarization Raman li- 
dar at 67.9°N in January 1997 reveal a strong positive correlation between the particle 
optical properties, specifically depolarization ratio <$ par and extinction-to-backscatter (li- 
dar) ratio S, for <S par <~ 40%, and an anti-correlation for <5 par >~ 40%. Over the length 
of the measurements the particle properties vary systematically. Initially, <5 par » 60% 
and 5 « 10 sr are observed. Then, with decreasing <5 par , 5 first increases to ~ 27 sr (<5 par « 
40%) before decreasing to values around 10 sr again (<S pa r ~ 20%). The analysis of li- 
dar humidity and radiosonde temperature data shows that the measured optical prop- 
erties stem from scattering by dry solid ice particles, while scattering by supercooled droplets, 
or by wetted or subliming ice particles can be excluded. For the microphysical interpre- 
tation of the lidar measurements, ray-tracing computations of particle scattering prop- 
erties have been used. The comparison with the theoretical data suggests that the ob- 
served cirrus data can be interpreted in terms of size, shape, and, under the assump- 
tion that the lidar measurements of consecutive cloud segments can be mapped on the 
temporal development of a single cloud parcel moving along its trajectory, growth of the 
cirrus particles: Near the cloud top in the early stage of cirrus development, light scat- 
tering by nearly isometric particles that have the optical characteristics of hexagonal columns 
(short, column-like particles) is dominant. Over time the ice particles grow, and as the 
cloud base height extends to lower altitudes characterized by warmer temperatures they 
become morphologically diverse. For large 5 and depolarization values of - 40%, the 
scattering contributions of column- and plate-like particles are roughly the same. In the 
lower ranges of the cirrus clouds, light scattering is predominantly by plate-like ice par- 
ticles. This interpretation assumes random orientation of the cirrus particles. Simulations 
with a simple model suggest, however, that the positive correlation between S and <5 par , 
which is observed for depolarization ratios < 40% mainly at low cloud altitudes, can 
be alternatively explained by horizontal alignment of a fraction of the cirrus particle pop- 
ulation. 


1. Introduction 

The microphysics of aerosols and clouds affect atmo- 
spheric processes of all scales [Baker, 1997]. The retrieval 
of microphysical information from remotely sensed optical 
data such as those acquired with lidar is therefore a worth- 
while research effort. If the sizes of the scattering parti- 
cles are comparable to the wavelength of the radiation used 
for the observation, one can exploit the wavelength depen- 
dence of the cloud optical properties to gain insight into the 
microphysics of the particulate scatterers. Over the years, 
sophisticated lidar instrumentation and retrieval algorithms 
have been developed for this purpose that have been suc- 
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cessfully applied to stratospheric aerosols [ Del Guasta et al. y 
1994; Wandinger et al, 1995], polar stratospheric clouds 
[Car slaw et al , 1998; Mehrtens et oh, 1999], tropospheric 
aerosols [Muller et al ., 1998; Muller et al., 2001], and con- 
trails [Sassen et al ., 2001]. 

However, these techniques are generally not applicable to 
tropospheric or stratospheric clouds with populations of par- 
ticles that are predominantly much larger than the optical 
observation wavelengths. For these clouds, other methods 
are required for the retrieval of microphysical information. 
Multi-instrument techniques, such as the approach based 
on combined lidar and radar measurements [Intrieri et al . , 
1993], have been proposed, but, since instrument clusters are 
not always available, it seems desirable to attempt extrac- 
tion of microphysical data from lidar measurements alone. 

For cirrus clouds, wdiich are the focus of this work, re- 
mote measurements of particle depolarization ratio (S par have 
proven to be highly valuable for studies of ice-cloud micro- 
physics, dynamics, and optical displays (see, e.g., Sassen 
[1991], Sassen [1992], Sassen and Takano [2000], Sassen and 
Benson [2001], Sassen et al [2001]). The success of using de- 
polarization ratio for the interpretation of cirrus data stems 
from the fact that <S par is one of the two cloud optical prop- 
erties accessible to lidar that do not depend on the particle 
number density. Thus, <5 pa r contains (particle population- 
averaged) information about particle morphology. To dis- 
tinguish <S par from cloud optical properties that do depend 
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on number density (such as extinction and backscatter co- 
efficient), we use the term particle optical property. 

The second particle optical property that can be mea- 
sured with lidar is the extinction-to-backscatter (lidar) ratio 
5, and it appears promising to include S in the microphysi- 
cal analysis of cirrus measurements. To our knowledge, only 
two cirrus studies based on lidar ratio and depolarization ra- 
tio have been published so far. Del Guasta [2001] retrieved 
cloud-averaged lidar ratios for Antarctic cirrus clouds from 
their measurements taken with an elastic polarization lidar, 
and interpreted the optical data in terms of ice crystal shape. 
Eloranta et al [2001] investigated correlations among parti- 
cle lidar and depolarization ratios of cirrus clouds measured 
over the continental United States with their polarization 
high spectral resolution lidar [ Eloranta and Pnronen , 1997], 
and compared experimental and theoretical values of lidar 
ratio. 

Recently, Reichardt et al. [2002b] have reported on the 
correlation between depolarization and lidar ratio found in 
Arctic and mid-latitude cirrus clouds, and discussed the im- 
plications of their findings for spaceborne remote sensing. In 
this follow-up article, we give a microphysical interpretation 
of these measurements. Since for the microphysical analysis 
it is instructive to look at the temporal evolution of the par- 
ticle optical properties of the cirrus clouds, we thus restrict 
our discussion to the long-duration cirrus observations made 
with the GKSS Raman lidar [Reichardt et al . , 1996] above 
the Swedish research facility Esrange (67.9 °N, 21.1 °E) in 
January 1997. 

2. Observations 

Figures 1-3 show three long-duration observations of cir- 
rus clouds made over Esrange in January 1997. These mea- 
surement cases have been chosen for outstanding data qual- 
ity; for a time-height visualization of the cirrus observations 
and a detailed description of the data evaluation process see 
Reichardt et al. [2002b]. For each lidar observation, three 
measurement intervals have been selected to study the de- 
velopment of the clouds over time. All cirrus clouds are op- 
tically thin, as can be seen from the particle and molecular 
backscatter coefficients. Generally, the particle backscatter 
coefficient exhibits a much larger variability with height than 
the particle depolarization ratio or the lidar ratio, indicat- 
ing that the particle number density rather than the particle 
optical properties, and hence the microphysical properties, 
changes rapidly within the clouds. 

In Fig. 4 the temporal evolution of the relation between 
the particle optical properties, specifically the particle de- 
polarization ratio and lidar ratio, is illustrated. The ex- 
perimental results can be summarized as follows. Firstly, 
for each measurement interval, at least one of the corre- 
lation patterns emerges that are evident in the combined 
data set of the three observation nights (Fig. 5). Secondly, 
during each cirrus event, the particle properties exhibit the 
same systematic variation over the length of each measure- 
ment (3-7.5 hours). Initially, J pa r ~ 60% and S ~ 10 sr 
are observed. Then, with decreasing 5 par , S first increases 
to ~ 27 sr (<5 par » 40%) before decreasing to values around 
10 sr again (5 par ~ 20%). 

Our data have been measured at a wavelength of 355 nm. 
Cirrus observations with the same instrument over north- 
ern Germany have yielded similar ranges of 5 and 5 par val- 
ues, except that <5 par > 35% is rarely measured [Reichardt, 
1999]. In contrast, if compared to other lidar data sets that 


Figure 1. Particle backscatter coefficient (/J), particle de- 
polarization ratio (^par)> and lidar ratio (5) for a cirrus sys- 
tem evolving over time. The measurement was taken over 
Esrange, Sweden, on 8 January 1997. Of the long-duration 
observation three measurement periods are shown (times 
corresponding to the centers of the intervals are indicated). 
Measurement wavelength is 355 nm. Error bars indicate the 
standard deviation due to signal noise. Molecular backscat- 
ter coefficients (dashed curves) are given for comparison. 


have been used for the retrieval of cirrus microphysics, some 
remarkable discrepancies can be found. Sassen and Corn- 
stock [2001] retrieved cloud- averaged lidar ratios at 694 nm 
for various categories of cirrus clouds observed over cen- 
tral North America. These lidar ratios tend to be larger 
than those shown in Figs. 1-4, the S distributions peak 
at 30 sr, and are weighted more heavily to larger S values. 
Likewise, Del Guasta et al. [1993] reported a considerable 
number of Antarctic cirrus cloud observations with cloud- 
averaged S > 30 sr for T c < -30°C (observation wavelength 
of 532 nm). Direct measurements of lidar-ratio profiles at 
532 nm by Eloranta et ai [2001] yielded a distribution of 
5 occurrences with greater similarity to our measurements 
(most of the data between 10 and 30 sr), but the center of 
the distribution appears to be at slightly larger S values 
than at 355 nm, and S data > 30 sr are observed (central 
North America site). On the other hand, particle depolar- 
ization ratios of Sassen et al. [2000] (£ par between 35% and 
60% for Tc < -25° C; tropical western Pacific Ocean; mea- 
surement wavelength of 1064 nm), Sassen and Benson [2001] 
(> 30% for T c < — 30°C) and of Eloranta et al. [2001] (most 
£ P ar > 30%) agree well with our observations, while the de- 
polarization data of Del Guasta et al. [1993] are significantly 
smaller with values predominantly < 22%. The discrepan- 
cies between the cirrus optical data sets mentioned may be 
due to the different geographical locations at which the cir- 
rus clouds were observed, and, in the case of lidar ratio, 
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Figure 2. Same as in Fig. 1, but for 16 January 1997. 


to the differing direct or indirect retrieval techniques. Dif- 
ferences in measurement wavelength might play a role also, 
although a wavelength effect would not be expected unless 
significant numbers of very small ice particles are present in 
the cirrus clouds. In view of this data review, it would be 
interesting to compare our observations with measurements 
obtained with other UV polarization Raman lidars at differ- 
ent locations. To our knowledge, however, no 355-nm data 
are available at this time. To address the issues discussed 
above, simultaneous measurements of depolarization and li- 
dar ratio at ultraviolet and visible wavelengths should be 
performed in the future. 

Figure 5 concisely summarizes the results of the corre- 
lation studies and the analysis of cirrus temperature asso- 
ciated with the particle optical properties which have been 
discussed by Reichardt et al. [2002b], and which are impor- 
tant for the microphysical interpretation proposed here. The 
combined data show two distinct linear correlation patterns. 
£par and 5 are negatively correlated when the depolarization 
ratio is high, and they are positively correlated when the de- 
polarization ratio is low. Furthermore, the particle optical 
properties depend on the ambient cloud temperature Tc* 
The coldest Tc values are found for high <5 par and low S, 
warmest Tc values for low <5 par and low S. For large 5, 
where the low- and high-J par data subsets merge, Tc lies 
approximately in the middle of the observed temperature 
range. Interestingly, this temperature is roughly the same as 
the mid-cloud temperature at which Platt and Dilley [1981] 
noted a change in their cirrus backscatter-to-extinction ra- 
tio data, a change they attributed to a transition in particle 
morphology. 

The correlations among S and <5 pa r, and their dependence 
on Tc suggest that the observations can be interpreted in 
terms of the shape and size of the scattering cirrus parti- 
cles (scatterers of any shape formed from frozen water) for 


Figure 3. Same as in Fig. 1, but for 19 January 1997. 


two reasons. Firstly, lidar-relevant particle optical proper- 
ties are determined by the phase matrix elements at 180- 
degree scattering angle of light scattering by the cirrus par- 
ticle population, which in turn depend on the morphology of 
the cloud particles. Secondly, Tc is the principal factor that 
controls basic crystal shape [Mason et al ., 1963]. Before we 
investigate this hypothesis, we first discuss the ray- tracing 
calculations we have performed in order to derive theoretical 
particle optical properties for comparison with the observa- 
tions. 

3. Ray-tracing Calculations 

Atmospheric crystals rarely have perfect shape, but ex- 
hibit distortions such as ice core inhomogeneities and angles 
between crystal faces that deviate from the theoretical val- 
ues. While particle depolarization ratios calculated with 
the ray-tracing technique often match values observed in 
cirrus clouds regardless of the assumed imperfections, the- 
oretical lidar ratios are very sensitive to the assumptions 
made and can be wrong by orders of magnitude. In order to 
derive realistic theoretical optical data, we account for crys- 
tal irregularities in our ray-tracing calculations by statistical 
roughening of the crystal surfaces [Hess et al ., 1998]. The 
crystals have principally hexagonal shape and are random 
in orientation, surface roughness is assumed to be the same 
for columns or plates of all aspect ratios. To reduce sta- 
tistical errors of the computations, scattering phase matrix 
elements calculated at 179- and 180-degree scattering angle 
are averaged. Model runs are performed for different sur- 
face roughnesses, and values of this parameter are selected 
for which the calculated <£ par and S values both fall within 
the observation ranges of these quantities. Figure 6 shows 
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Figure 4. Particle depolarization ratio versus lidar ratio 
scatter plots of the cirrus measurements on 8 (top), 16 (cen- 
ter), and 19 (bottom) January 1997 shown in Figs. 1-3. 


the results of calculations for hexagons with 0%, 1%, and 
10% surface roughness. Generally, an increase in surface 
roughness leads to higher S and lower <5 pa r- For the compar- 
ison with our observations, we choose data obtained with 
0% (plates) and 1% (columns) surface roughness. 

One shortcoming of the ray-tracing approach is that 
physical-optics effects are not taken into account. Fur- 
thermore, numerical instability is a problem in the direct 
backward-scattering direction, so ray-tracing data have to 
be used with caution if lidar measurements are to be inter- 
preted. For crystals with imperfections, however, the ray- 
tracing technique may still be applicable, because crystal 
shape distortions tend to smooth the features of the scatter- 


Figure 5. Combined particle depolarization ratio versus 
lidar ratio scatter plot of the cirrus measurements shown 
in Figs. 1-3. The distribution of the data points is bifur- 
cated. The domains chosen for splitting the data set into 
two subsets of data points with high <5> par and low <$ par are 
indicated by differently slanted lines. Within each group, 
data are strongly correlated. Corresponding regression lines 
are shown. Cloud temperatures that are associated with the 
particle optical properties of each data subset are indicated 
for lidar ratios of 12 and 25 sr (cf. Reichardt et al. [2002b]). 


Figure 6. Particle depolarization ratio versus lidar ratio 
scatter plot of solid hexagons with aspect ratios between 
0.04 and 0.3 (plates; A, A), and aspect ratios between 1.3 
and 5.9 (columns; V, ▼ ). The data have been derived with 
ray-tracing calculations. Assumed surface roughness of the 
artificial crystals is indicated, filling of symbols alternates 
with surface roughness. The display range of the cirrus ob- 
servations is shown for comparison (dotted rectangle). 

ing phase matrix elements at lateral and backward scattering 
angles. 

4. Microphysical Interpretation 

Figure 7 shows the domains of the observed particle op- 
tical data along with the results of the ray-tracing calcula- 
tions. Depolarization ratios of relatively isometric crystals 
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Figure 7. Theoretical optical data of plates (A) and 
columns (▼), and domains and regression lines of the high- 
5 par and low-5 par data subsets. Aspect ratios of plates de- 
crease, aspect ratios of columns increase with lidar ratio. 


match the values of the high-5 (cold 7c) and low-5 (warm 
7c) data subsets for small S. For larger S, the computed op- 
tical properties of elongated columns and thin plates mark 
the boundaries of the measurements. This implies that for 
large S the observed cirrus clouds have the optical proper- 
ties of mixtures of columnar and planar particles. Since in 
natural ice clouds length and aspect ratio (length divided by 
width) of the ice crystals are correlated (small crystals tend 
to have aspect ratios close to one, whereas large columns and 
plates have aspect ratios 1 and <C 1, respectively [Auer 
and Veal , 1970]), the optical data indicates an increase in 
size of the scattering cirrus crystals with 5. 

Despite the good agreement between theoretical and ob- 
served optical properties, however, the translation of the 
optical properties into actual shapes and sizes of the cirrus 
particles is not straight forward. It is complicated by the fact 
that cirrus clouds usually contain particles of various sizes 
and shapes, and that the relation between microphysical and 
optical properties is not unique. E g., hexagonal plates and 
hollow columns, bullet rosettes, and aggregates have almost 
the same scattering characteristics, if rough crystal surfaces 
are assumed [Yang and Liou, 1998]. For this reason, it would 
be more accurate to say that plate-Jt&e or column- like parti- 
cles dominated cirrus scattering when the optical signature 
of hexagonal plates or columns was observed. Furthermore, 
composite particles have been shown to scatter similarly to 
the ensemble of basic crystals that formed it [Macke, 1993; 
Yang and Lion , 1998]. So a small ice particle consisting 
of one or more relatively large ice crystals may still have 
a higher S value than a large ice particle formed from rela- 
tively small crystals. With these limitations in mind, we pro- 
pose the following microphysical interpretation of the obser- 
vations. In the coldest regions of the clouds, approximately 
isometric ice particles form. For increasing 7c, the ice par- 
ticles get larger and consist of a mixture of column-like and 
plate-like scatterers, but remain predominantly column-like. 
Around To = -50°C, the dominant characteristics of the 
particles gradually shifts from column- to plate-like. At 
even warmer temperatures in the lower regions of the cirrus 
clouds, size and morphological complexity of the plate-like 


particles increase with 7fc. This interpretation combines in 
a consistent manner the observed cirrus optical properties 
and temperature data with theoretical studies on light scat- 
tering by atmospheric ice particles [Macke, 1993; Hess et 
a/., 1998; Yang and Liou , 1998], laboratory experiments on 
the formation and growth of ice crystals [ Heymsfield , 1973; 
Gonda and Yamazaki , 1978], and in situ particle sampling 
[. Heymsfield , 1973]. According to the proposed scenario, the 
ice clouds observed over Esrange would have had a vertical 
habit distribution that followed the common, yet not gen- 
eral, pattern of small, sphere-like crystals near the cloud top, 
relatively large mono-crystals of imperfect hexagonal sym- 
metry or simple bullet rosettes around the center of the ice 
cloud, and crystal clusters of arbitrary over-all shape in the 
lower cirrus ranges (see Heymsfield and Iaquinta [2000] and 
Yang et al. [2001] for good examples of height-dependent 
habit variations in cold cirrus clouds). 

It is important to note that this microphysical interpre- 
tation does not draw upon the observed evolution of the 
optical properties over time, as shown in Fig. 4. It is a 
static interpretation in the sense that for each of the three 
cirrus cases considered all observation intervals are treated 
as if independent of one another. This assumption led to 
the time-independent presentation of the experimental data 
in Fig. 5, and to the interpretation given above. However, 
Fig. 4 strongly suggests that the systematic change in par- 
ticle optical properties from one measurement interval to 
the other may be interpreted in terms of growth; or more 
generally morphological development of the cirrus particles. 
The difficulty with this dynamical interpretation of the lidar 
data is that one has to assume that the lidar observations 
of consecutive cloud segments can be mapped on the tem- 
poral development of a single cloud parcel moving along its 
trajectory. This requires the cloud layer, though varying in 
vertical structure and particle number density, to be com- 
posed of similar particles throughout its horizontal range 
at any given time over the observation period. Because of 
the generally high variability of the atmosphere, it cannot 
be expected that this condition is often met, and therefore 
the dynamical interpretation is a more speculative one than 
the static interpretation at this point. However, at least in 
one measurement case, initial results obtained with a cou- 
pled microphysical and optical model currently being devel- 
oped are consistent with the hypothesis of particle growth 
[Reichardt et al 2002a]. 

5. Discussion 

The interpretation given above is plausible. However, 
since extraction of microphysical information from optical 
remote measurements is an ill-posed problem, and since at- 
mospheric ice particles can have diverse morphology, alter- 
native explanations are conceivable. E.g., the decrease in 
particle depolarization ratio with S observed in the high- 
5 P ar data subset may not be due to an increasing fraction 
of plate-like particles. Instead it may be the effect of col- 
umn basal faces hollowing out during growth, since hollow 
columns have lower 5 par values than solid specimen [Takano 
and Liou , 1995]. We have tried to avoid this morphologi- 
cal problem by referring to the cirrus particles as column- 
or plate-like particles rather than columns or plates. We 
may be able to refine our microphysical interpretation when 
reliable theoretical optical data for ice particles other than 
columns and plates become available to us. Here, except for 
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shape, we discuss the other underlying fundamental assump- 
tions of our microphysical interpretation. Specifically, these 
assumptions are that the cirrus particles consist of frozen 
water and are in random orientation, and that sublimation 
does not play a role. 

5.1. Cirrus Particle Phase 

Because of the very cold cloud temperatures (below 
— 40° C) for the overwhelming majority of the cirrus data), 
it is a justified assumption that scattering is by dry solid 
ice particles only. At these temperatures, noticeable effects 
are highly unlikely that stem from quasi-liquid layers on the 
crystals [Wettlaufer, 1999] or that are due to co-existing liq- 
uid droplets. To our knowledge, the coldest temperatures at 
which liquid droplets have been reported to exist in cirrus 
and convective clouds are -36° C and -37.5° C, respectively 
[Sassen et a/., 1985; Rosenfeld and Woodley , 2000]. Liquid 
droplets of aqueous acids can exist at even colder atmo- 
spheric temperatures. Freezing point depression resulting 
from an influx of stratospheric volcanic aerosols into the 
cirrus-generating regions has been used to explain the ob- 
servation of short-lived cloud cells of liquid drops in cirrus 
clouds at -49° C [ Sassen , 1992]. However, this rare phe- 
nomenon can be excluded from the discussion here, because 
at the time of the measurements the aerosol content of the 
lower stratosphere was very low. 

5.2. Random Orientation 

The microphysical interpretation given in Sect. 4 is based 
on the assumption that at all measurement times the cirrus 
particles were randomly oriented in three dimensions (3D- 
orientation). As a consequence, we had to invoke an increase 
in morphological complexity of plate-like particles in order 
to explain the low 5 values in the lower cirrus layers. 

An alternative mechanism that could explain more ele- 
gantly the low-^par data subset, as well as (if the dynamical 
interpretation is considered) the transition between the two 
groups of correlated data, is the departure of an increasing 
fraction of the particle population from a random orienta- 
tion as one approaches the cloud base. Falling ice parti- 
cles larger than a critical size (> 100-200 /an, depending on 
shape [Sassen, 1980]) have the tendency to assume an ori- 
entation with their maximum geometrical dimension in the 
horizontal plane (2D-orientation). If crystal faces of these 
2D-oriented ice particles are (almost) perfectly aligned hori- 
zontally, the lidar signature shows extremely low <$ pa r and S 
values due to specular reflections of laser light into the lidar 
receiver. However, it is relatively unlikely that 2D-oriented 
columnar particles are aligned in this peculiar manner, since 
they can rotate about their horizontally oriented long axis 
easily. This is the reason why the optical display associ- 
ated with this special orientation of hexagonal columns, the 
Parry arc, is rarely observed [Sassen and Takano , 2000]. As 
a consequence, horizontal alignment of columns remains gen- 
erally undetected by lidar. Contrarily to the prism faces of 
aligned columns, horizontal orientation of the basal crystal 
faces of 2D-oriented plates is better preserved, even though 
they may fall with a swaying motion. Therefore it is gen- 
erally accepted that observed very low depolarization ratios 
[Platt et ai , 1978] and lidar ratios [Ansmann et al, 1992] 
are generated by 2D-oriented planar particles [Platt, 1978; 
Sassen , 1980; Sassen and Benson , 2001]. 

In the following, we show that the assumption of 2D- 
oriented ice particles being present is consistent with the lin- 
ear correlation between S and (5 par observed in the low-5 par 
cirrus data. In our simple model, we assume that the sim- 
ulated cirrus cloud consists of a population of identical ice 


Figure 8. Profiles of particle backscatter coefficient, parti- 
cle depolarization ratio, and lidar ratio for the cirrus cloud 
measured over Geesthacht, Germany, around 17:23 UT on 
27 November 1995 (top row). Values of 0 below 7 km have 
been multiplied by 0.33 to fit into the display range. Parti- 
cle depolarization ratio versus lidar ratio scatter plot of this 
measurement, and calculated optical data of a hypothetical 
cirrus cloud consisting of randomly oriented (5 = 29.5 sr, 
<5 par = 40%) and horizontally aligned (5 = 2sr, <5 par — 2%) 
ice particles (thick solid curve) (bottom). Numbers indicate 
fraction of aligned particles. Display range of the cirrus 
observations over Esrange, and data-subset domains and re- 
gression lines are shown for comparison. 

The measurement on 27 November 1995 is one of the rare 
cirrus cases for which reliable height-resolved profiles of lidar 
and depolarization ratio could be obtained over Geesthacht 
(because measurement conditions were generally not favor- 
able). <5 P ar and s 3X6 correlated; for S > 8sr, data points 
fall into the Esrange low-depolarization domain. Predomi- 
nantly, these data points come from the cloud layer above 
8 km. 


particles, of which a fraction a is in perfect 2D-orientation 
while the other particles remain randomly oriented. Then 
the cloud optical properties as a function of a are given by 


5(a) - 

^par(O') ~ 


1 + q(i^G ~ 1) £3D 

l+a(R G Rs l -l) 

i+a{^i(i+0(i+C)W -n f 3 D 

1 + a{[(l + *$)/(! + 6 2 P Z)}R G Rs l - U ^ 


where Rq = G 2D /G 3D , R$ = S 2D /S 3D , and R$ — <5p? r /^par 
are the ratios of geometrical cross section G , lidar ratio, and 
depolarization ratio, respectively, for 2D-oriented cirrus par- 
ticles with respect to the randomly oriented 3D case. For 
large thin plates, we have R G ~ 2. This means that on 
horizontal orientation, G of the particles, and hence the ex- 
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tinction coefficient, increase by this relatively small factor. 
On the other hand, the extinction-to-backscatter (lidar) ra- 
tios of 2D-oriented particles are 5-10 times smaller than 
those of randomly oriented scatterers. Apparently, even in 
cirrus clouds with small fractions of 2D-oriented particles, 
cloud backscattering is strongly affected by the specular re- 
flections of the aligned particles, while cloud extinction is 
still determined by the entire particle population. 

In order to calculate S(a) and <S pa r(a) 3 we nee< ^ to know 
reliable values for lidar and depolarization ratios of cirrus 
particles in 2D- and 3D-orientation. Optical modeling can 
provide S 3D and 5™ , however, due to methodological lim- 
itations, calculations for 2D-oriented particles are currently 
not available. For this reason, we use instead as input to 
our calculations experimental particle optical data that can 
be clearly attributed to light scattering by ice particles in 
either random or horizontal orientation. For 3D-oriented 
particles we assume S 3D = 29.5 sr and (S p ? r = 40%, which 
are approximately the optical properties of particles in the 
low- <5 par data group for large 5. 

Values for S 2D and 8™ r are taken from cirrus data ac- 
quired on 27 January 1995 over Geesthacht, Germany, with 
the same lidar instrument that was deployed in Esrange. 
Cloud profiles of backscatter coefficient, particle depolariza- 
tion ratio, and lidar ratio are shown in Fig. 8. On this day, 
two layers of ice clouds were observed. In the lower layer, 
J par and 5 decreased gradually from the cloud base to the 
height of the backscatter maximum, and increased subse- 
quently to the cloud-layer top. <5 pa r and 5 at the boundaries 
of the lower cloud layer are similar to those found in the 
upper cloud layer. This optical signature is certainly the 
result of horizontal particle orientation, and we choose the 
minimum values of lidar and depolarization ratio, i. e., 2sr 
and 2%, respectively, for S 2D and J P a r in our simulations. 

The model results are presented in Fig. 8, along with the 
domains of the Esrange data subsets. As expected, align- 
ment of a relatively small fraction of cirrus particles leads 
to a substantial reduction in S and <5 par - The synthetic de- 
polarization ratio exhibits a dependence on lidar ratio that 
is similar to the measured correlation. Considering the sim- 
plicity of the model, it is not surprising that simulated and 
observed (S pa r- versus - S data do not match exactly. With 
<5 p ar ~ 10%, lt is P ossible t0 achieve close agreement be- 
tween model and measurement, but sensitivity studies show 
that the theoretical relation only slightly deviates from the 
one shown, if S 2D and <5 P ar are varied within limits that are 
reasonable for aligned planar cirrus particles 2-5 sr, ~ 2- 
5%). This study supports the idea that horizontal particle 
alignment is the main cause for the relation between S and 
<$ P ar in the low-depolarization domain of the Esrange data. 
However, as detailed below, there are arguments against this 
interpretation. 

If orientation of the cirrus particles occurred during our 
measurements over Esrange, why then did we never detect 
extremely low lidar and depolarization ratios such as those 
observed in the lower cirrus layer of the Geesthacht measure- 
ment? One explanation could be that the GKSS Raman 
lidar was not sensitive to light scattering by horizontally 
aligned cirrus particles at that time. This mobile lidar sys- 
tem was transported to Esrange in December 1996, and was 
setup for its deployment under difficult winter conditions. 
So it is perfectly possible that the instrument was point- 
ing a couple of degrees off zenith. In this case, as shown 
by Sassen and Benson [2001], the influence of 2D-oriented 
ice particles on the lidar measurements would have been 


Figure 9. Profiles of particle backscatter coefficient, parti- 
cle depolarization ratio, and lidar ratio for the cirrus cloud 
measured over Esrange, Sweden, around 05:04 UT on 15 
January 1997 (top row), and particle depolarization ratio 
versus lidar ratio scatter plot of this measurement (bottom 
row). Theoretical optical data, data-subset domains and 
regression lines are shown for comparison. 


greatly reduced, if not annihilated, and the Esrange data dis- 
cussed here would have the optical signature of (quasi) ran- 
domly oriented scatterers. Another reason for the absence of 
S < 8sr and 6 pAT < 15% in the Esrange data could be that 
the oriented particles sensed by the lidar were in themselves 
disordered aggregations of crystals. To our knowledge, scat- 
tering signatures of these particles have not been calculated, 
but it can be speculated that the optical properties of these 
particles, if aligned in the horizontal plane, would be similar 
to those of the randomly arrayed ensemble of constituting 
crystals, thus covering up their reduced orientational free- 
dom [Sassen, 1980]. If this were true, we would have to 
infer that the 2D-oriented cirrus particles which produced 
the optical display of extreme particle properties over north- 
ern Germany had a simpler morphology than those observed 
above Esrange. In view of the growth conditions at Esrange 
(from lidar humidity and radiosonde temperature data the 
diffusional mass growth rate [Pruppacher and Klett, 1980] 
can be estimated to be at least a factor of two higher over 
Geesthacht than over Esrange), this appears to be rather 
unlikely. 

With regard to the issue of particle orientation, the dis- 
cussion of the cirrus measurement taken on 15 January 1997 
over Esrange is instructive (Fig. 9). On this day, depolar- 
ization ratios decreased with height without much variation, 
with (S P ar values between 24% at the base and 16% at the top 
of the cloud. The observation of high values of S (20-27 sr) 
makes it impossible to explain these rather low depolariza- 
tion ratios in terms of particle horizontal alignment. Instead, 
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the agreement between measured and theoretical data sug- 
gests that the ice cloud on 15 January 1997 consisted of 
large plate-like particles in 3D-orientation. This cirrus mea- 
surement supports the underlying assumption about orien- 
tational randomness that is made for the microphysical in- 
terpretation put forward in Section 4. However, it also shows 
that not in all measurement cases the optical properties of 
cirrus clouds with large lidar ratios carry the signature of 
mixtures of plate-like and column-like ice particles (on 15 
January 1997, S pAT and 5 were correlated, but column-like 
scatterers were apparently absent). 

In summary, the correlation between lidar and depolar- 
ization ratio found for S pAT <40% may be viewed as the 
result of horizontal alignment of planar ice particles, or may 
be attributed to the presence of randomly oriented plate-like 
scatterers. Since good arguments exist for the two positions 
(or for a combination of both), more studies are required to 
shed some light on the mechanisms that lead to this optical 
phenomenon. 

5.3. Sublimation 

It might be argued that sublimation led to the observed 
decrease in depolarization ratio for warmer temperatures 
and lower cloud altitudes. Experimental studies show that 
evaporating ice crystals acquire rounded shapes with smooth 
surfaces [ Nelson , 1998], and it appears obvious to con- 
clude that these increasingly spheroidal particles change 
their scattering behavior toward lower depolarization ra- 
tios. However, this argument is not conclusive. Scattering 
phase matrices are predominantly determined by reflections 
at crystal surfaces, and therefore sublimation-induced alter- 
ation of edge and corner radii may be better regarded as dis- 
tortions of the initial particle habits, having relatively little 
effect on the scattering signature, specifically the depolariza- 
tion ratio. Laboratory and theoretical studies seem to sup- 
port this assumption. Sassen [1977] reported depolarization 
ratios between 52% and 59% for frozen rainwater drops with 
a regular spherical or spheroidal appearance. T-matrix cal- 
culations performed for spheroids, circular cylinders, and de- 
formed spheroids (Chebyshev particles) by Mishchenko and 
Sassen [1998] show that these rounded crystals have depo- 
larization ratios that are comparable with, or larger than, 
<$ P ar values of hexagonal plates obtained with the ray-tracing 
method. Although the size parameters of the particles as- 
sumed for the T-matrix calculations are smaller than those 
normally found in cirrus clouds, these theoretical results, 
just like the experimental findings, suggest that evapora- 
tion of the ice crystals does not have to be invoked for the 
explanation of the observed low depolarization ratios. Sub- 
limation may be also ruled out as a major cause for the low- 
depolarization data subset for observational reasons. Simul- 
taneously with the cirrus measurements, water vapor was 
measured with the GKSS Raman lidar. From the moisture 
data and from radiosonde temperature observations we re- 
trieved estimates of relative humidity with respect to ice 
(Fig. 10). The low absolute humidity over Esrange resulted 
in weak water-vapor lidar signals that require long data in- 
tegration times and restrict this analysis to altitudes < 6- 
7 km. This does not pose a major problem for the discussion 
here, because measurement segments with relatively stable 
lower cloud boundary could be chosen for the integration, 
and most of the low-J par cirrus data are found towards the 
bases of the clouds. A more serious limitation of the quality 


Figure 10. Relative humidity with respect to ice (RHI) re- 
trieved from lidar water-vapor and radiosonde temperature 
data. Lidar data were integrated over measurement inter- 
vals for which the cirrus clouds extended to low altitudes 
(times corresponding to the centers of the intervals are indi- 
cated). Range of cloud base height during the water- vapor 
measurement intervals (shaded areas), and cloud base height 
as observed during the cirrus measurement segments with 
low <5par (horizontal solid line) are given for comparison. Er- 
ror bars indicate the standard deviation due to signal noise. 
Temperature profiles were taken from radiosondes measure- 
ments at the lidar site (8 January 1997, 14:00 UT), or at the 
radiosonde station Lulea (240 km south of Esrange; 17 and 
20 January 1997, 00:00 UT). 


of the relative-humidity data that complicates its interpre- 
tation is that temperature was not observed together with 
the lidar measurements at the lidar site. Although caution 
has to be applied, the relative-humidity data indicate that 
in all cases the cirrus clouds developed in air layers which 
were most likely saturated or supersaturated with respect 
to ice (0-20% supersaturation) at cirrus altitudes where low 
particle depolarization is found. This conclusion is corrobo- 
rated by the analysis of the cloud extinction profiles, as the 
extinction coefficient is usually at its maximum in the lower 
cirrus layers, indicating particle growth rather than evapo- 
ration. For these reasons, sublimation is not considered in 
our microphysical interpretation. 

6. Conclusions 

We have shown in this study how the polarization Raman 
lidar data of cirrus clouds that evolved over the Swedish re- 
search facility in January 1997 can be interpreted in terms 
of the sizes and shapes of the scattering ice particles, and, 
under additional assumptions about the atmospheric condi- 
tions prevalent during the observations, in terms of growth, 
or more general, morphological development of the cirrus 
particles. It turned out that key to the understanding of 
light scattering by cirrus clouds is simultaneous measure- 
ment of particle depolarization ratios and lidar ratios. How- 
ever, we would like to stress the obvious fact that, although 
the given interpretation appears plausible, it is by no means 
proven. Important issues remain to be resolved which in- 
clude, among others, the relation between the optical and 
the microphysical properties of realistic ice particles, the in- 
fluence of particle orientation, and the dynamic interpreta- 
tion of the data. Furthermore, the question of the generality 
of our findings has to be addressed. 

In order to improve the microphysical analysis of cirrus 
lidar data, more research, both experimental and theoret- 
ical, has to be conducted. Two especially helpful experi- 
ments would consist of concurrent lidar observations from 
the ground and airborne particle sampling, and simultane- 
ous in situ measurements with microphysical instrumenta- 
tion and a lidar analog. We have demonstrated that knowl- 
edge of supersaturation aids the interpretation of the cirrus 
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data. A shortcoming of our supersaturation data was that 
temperatures were not measured with the lidar, but had to 
be taken from radiosonde observations. Addition of cloud- 
temperature measurement capability to a lidar deployed in 
cirrus observations would be certainly desirable; the method 
that can be used for cloud temperature studies is the rota- 
tional Raman technique [Behrendt and Reichardt , 2000]. 

Also, improvements in calculation of optical properties 
of large crystals are necessary. The ray-by-ray-integration 
method appears to be well-suited for this task, since some of 
the restrictions that standard ray-tracing codes have do not 
apply here [Yang and Lion , 1997]. Finally, coupled optical 
and microphysical modeling might contribute to the under- 
standing of the microphysical information content that is 
hidden in lidar measurements of cirrus clouds [Reichardt et 
a/., 2002a], 
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